Theory of phonon-assisted "forbidden" optical transitions in spin-gapped systems 
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We consider the absorption of light with emission of one Stot = 1 magnetic excitation in systems 
with a spin gap induced by quantum fluctuations. We argue that an electric dipole transition is 
allowed on the condition that a virtual phonon instantaneously breaks the inversion symmetry. 
We derive an effective operator for the transition and argue that the proposed theory explains the 
polarized experiments in CuGeOa and SrCu2 (603)2 . 



I. INTRODUCTION 



Techniques of using the interactions between hght and 
spin-waves to study the excitations of magnetic sohds 
were developed shortly after the invention of the laser. 
Single magnon scattering of photons was first predicted 
from the Zeeman coupling of the magnetic field of the 
photon field to the magnetic spins, leading to mag- 
netic dipole transitions.^ Later it was pointed out^ 
that the electric field of the electromagnetic radiation 
could also couple to the spin, by an indirect process 
in which spin-orbit interactions act on electronic states 
excited virtually by electric-dipole transitions. Experi- 
ments in antiferromagnetsA. showed that this latter mech- 
anism dominated the magnetic-dipole transitions to sin- 
gle magnon excitations. The Raman spectrum also re- 
vealed relatively strong two-magnon scattering. This was 
argued^ to be due to an independent mechanism: excited- 
state exchange interactions. The same mechanism, by 
which the magnetic exchange interaction is modified by 
electric-dipole excitation of the magnetic electrons, was 
advanced^ to explain far-infrared absorption. A variant 
is to replace the virtual electronic excitation by a virtual 
lattice distortion that modifies the magnetic exchangei^ 
The intensities of such transitions can be calculated by 
writing effective operators for absorption or Raman scat- 
tering in terms of the spin operators^ii This theory is 
considered generally to give good account of inelastic 
light scattering and optical absorption. For an isotropic 
system the effective operator conserves total spin and 
what is commonly called the "Fleury-Loudon" theory is 
used to analyse the spectroscopy of spin conserving tran- 
sitions. 

Optical techniques are now well established as probes 
of magnetic excitations, whether it be by Raman scatter- 
ing, i.e. inelastic scattering of optical frequencies, elec- 
tron spin resonance (ESR), i.e. resonant absorption of 
electromagnetic radiation with sweeping magnetic field, 
or by transmission measurements of infrared radiation. 
The techniques have been further enhanced by the in- 
creasing flexibility of light sources and detectors in the 
far-infrared region that is useful to much of magnetism. 
ESR studies using sources derived from far-infrared lasers 
rather than the traditional cavities are now available up 
to THz frequencies and may be made in large static 



or pulsed magnetic fieldsi^ Transmission studies in the 
far infrared range have the advantage of allowing for 
measurement in zero external magnetic field. While re- 
stricted to small momentum transfer, (7 w 0, compared 
to neutron inelastic scattering, the optical techniques 
have the advantage of much higher frequency resolution. 
The possibility of polarising the electromagnetic radia- 
tion means different transition mechanisms may be dis- 
tinguished. 

Optical measurements are particularly useful for pre- 
cise measurements of the spin gap properties in strongly 
correlated systems and spin-liquid systems with mag- 
netic singlet ground states. Because of the frequencies 
now available, one can apply an electromagnetic source 
with sufficient energy to excite the first triplet Stot = 1 
excited state from the singlet Stot = ground state. 
Many systems of interest are highly isotropic with re- 
spect to spin rotations and transitions between the sin- 
glet Stot = ground state of the spin-liquid to the 
first triplet Stot — 1 excited state would be forbidden 
by symmetry in the isotropic limit. Even the weaker 
magnetic-dipole coupling should give zero intensity as 
the ground state is a spin singlet. One would then ex- 
pect to see the excited singlets, i.e. two magnon states 
only. Nonetheless the "forbidden" transitions to the sin- 
gle magnon states have been observed in many spin- 
liquid, ranging from the S=l/2 quasi-one-dimensional 
systems CuGeO^SiiaiiiiiSiiiii and NaVzOgi^, to 2d 
system such as SrCu2(B03)242ii^ and to the spin-1 chain 
compound, NENP^^. Despite detailed experiments, no 
clear understanding of the mechanism of these transi- 
tions has emerged. It is clear that spin-orbit coupling, 
that breaks the conservation of total spin, must be in- 
cluded as it is then possible a priori to have a transi- 
tion to a one-magnon state. As mentioned, the photon 
can couple to the spin degrees of freedom in different 
ways, via direct magnetic dipole transitions or indirect 
electric dipole transitions with spin-phonon or spin-orbit 
couplings. As one of the purposes of performing high res- 
olution spectroscopy is to resolve the weak anisotropics, 
it is important to distinguish between these mechanisms, 
i.e. to find the one which gives the strongest absorption. 
As in the original studies'* this is done by establishing, 
and then verifying experimentally, selection rules. For 
one-magnon absorption, previous estimations favored a 
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purely electric dipole transition for NENPr-- In the case 
of CuGeOa the suggestion that a staggered field would 
give rise to a magnetic dipole transition^ has been ruled 
out by the polarized experiments.^^ Furthermore the first 
order corrections to the Hamiltonian in spin-orbit cou- 
pling lead to vanishing magnetic dipole intensity owing to 
a lattice selection rule.^^ In the compound SrCu2 (603)2 
it has been shown experimentally that varying the direc- 
tion of the electric field of the wave (while keeping the 
magnetic field of the wave fixed) changes the intensity of 
the absorption, suggesting that the transition is electric- 
dipole in naturei—^ One would also like to know which 
of the two electric dipole mechanisms applies, absorp- 
tion involving solely the electronic degrees of freedom or 
with the lattice degrees of freedom. In the original the- 
ory of Elliott and Loudon of light scattering by magnons, 
the electric dipole coupling indeed leads to the creation 
of one-magnon excitations.^'* Although such two photon 
processes are not forbidden in infrared absorption, they 
are much smaller in intensity since they involve the weak 
coupling to light to second order in perturbation the- 
ory. Alternatively in the presence of strong spin-orbit 
coupling, it is possible to have single photon coupling 
to spin excitations^Si but as this is of second order in 
the spin-orbit coupling, we shall assume that the linear 
order will dominate for these materials, which are close 
to isotropic. In addition lattice symmetries such as cen- 
ters of inversion between the magnetic ion may eliminate 
such terms, or at least reduce them further, if the in- 
version symmetry is slightly broken, as in SrCu2 (603)2 
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In this paper we shall show that an effective operator of 
Dzyaloshinski-Moriya symmetrji^SiSi acting on the spin 
degrees of freedom, 

He^ E'3(i)A;37(a)(S. X S,+a)T (1) 

^,a,/3,7 

can be used to explain the polarized experiments of 
CuGe03 and SrCu2(B03)2 . Here E'^{t) is the com- 
ponent P of the applied electromagnetic field at time t. 
The indices i and a define the lattice of magnetic bonds 
and the coefhcients A./Sj will be made explicit in section 
im They couple the component P of the electric field with 
the component 7 of the vector product of the spin oper- 
ators. An electric dipole operator (1) can arise from an 
electronic mechanism, as may be the case in NENP— , but 
centers of inversion at the middle of the Cu-Cu bonds in 
CuGeOa and SrCu2(B03)2^'*, would forbid generation of 
the operator from purely electronic processes. A lattice 
distortion may, however, break the inversion symmetry 
instantaneously, and allow terms of the form in (1). We 
therefore consider the phonons explicitly, and in section 
nil we derive in detail the effective transition operator, 
which includes an anisotropic part of the form (1). The 
essential physical mechanism is that the electric field ex- 
cites a virtual phonon state Stot — which is coupled to 
the Stot = 1 state by an anisotropic spin-phonon coupling 



which originates in spin-orbit coupling. An explanation 
involving the modulation of static Dzyaloshinski-Moriya 
interactions has been put forward recently for the case 
of NaV205 ."^ In that compound, however, no polarized 
experiments are available and moreover, it is difficult to 
distinguish with a magnetic dipole transition which turns 
out not to be forbidden by a lattice selection rule.^^ The 
mechanism we develop here is more general in that it does 
not require the presence of a static Dzyaloshinski-Moriya 
interaction. It only needs the instantaneous breaking of 
the inversion center which is assured by the appropri- 
ate phonons. This allows us to consider the operator 
on the strongest bonds irrespective of whether the bond 
lacks an inversion center or not. In section^ we give the 
selection rules and the order of magnitude of such electric 
dipole transitions. We compare with the experiments in 
CuGeOa and SrCu2(B03)2 in section Iml 

II. EFFECTIVE MAGNETIC OPERATOR AND 
SELECTION RULES 

In this section we show that the first-order spin-orbit 
correction to the spin-phonon coupling leads indeed to 
an effective magnetic operator for the optical transitions. 
We note that a phonon- assisted optical transition is the 
usual explanation for the occurrence of the singlet Stot = 
bound states of two magnon states in the spectrum of 
the high Tc's cuprates.^ The spin-orbit correction should 
then lead to transition to Stot = 1 states. 

We start with a magnetic Hamiltonian for a chain or 
a layer of Cu atoms, for instance, that can be motivated 
by the usual super-exchange arguments: 

ff-^S,J({u,4)S,+a+i^ph-E.Pp;. (2) 

iad 

where Si is a spin operator, u^^ is the displacement vector 
of the ion d in the unit-cell i, Hph is the phonon Hamilto- 
nian which takes into account the kinetic part of the ions 
and the spring constants, Pph is the electric dipole of the 
ions and E is the external electric field. The magnetic 
couplings, 3{{uid}), can be expanded to ffist order in the 
ion displacements. Including the first order in spin-orbit 
coupling, there is an extra term of Dzyaloshinski-Moriya 
symmetry: 

Hsp = 9>?dS^■S.+a + dfu^d{S, X S,+,f (3) 

iada[5 

where is the partial derivative of the diagonal part of 
J({uid}) with respect to u^^ (it depends on the bond i, a 
but we will not write it explicitly in the following). The 
origin of d!^^ is explained below. This is indeed a general 
form for the spin-phonon coupling and there is no restric- 
tion to be added on the grounds of symmetry. The static 
Dzyaloshinski-Moriya interaction is forbidden when there 
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is an inversion center at the middle of the bond. If the set 
of displacements u^^^ is such as to remove the inversion 
center (which is the general case) then such an interaction 
takes place. For example if we take the two symmetric 
ninety degrees super-exchange paths Cu-O-Cu, there is a 
center of inversion and there is an interference between 
the two paths that leads to no Dzyaloshinski-Moriya in- 
teraction. Suppose now that the two oxygens move up- 
wards. Because the hopping of the electrons is much 
faster than the typical phonon frequency, the electrons 
see a frozen distorted lattice on that time scale. The 
interference therefore does not occur anymore and there 
is an effective Dzyaloshinski-Moriya interaction linear in 
the displacements in the first order. This is the origin of 
the second term of Q which involves a tensor since 
the displacements in one direction, a, generally produce 
a Dzyaloshinski-Moriya vector in another direction, f3. 
Strictly speaking, d^'^ also depends upon the bond i,a, 
but we do not write it explicitly. Note that this term 
is derived in a super-exchange approach by taking into 
account the spin-orbit coupling in first-order in perturba- 
tion theory in the lines of the original Moriya's articlei^ 
We shall refer to it as a dynamical Dzyaloshinski-Moriya 
interaction in the following. 

The transition probability is then given at zero tem- 
perature by the "golden rule" : 



Pph 



= I (/|E.Pp,|0) I'Siu;- 



(4) 
(5) 
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where w/ is the energy of the excitation, typically the 
one-magnon energy. At first order in Hgp in perturbation 
theory the matrix element is written in terms of a sum 
over the excited states: 



n 

The intermediate states that contribute to the sum over n 
contain one phonon (whereas the initial and final states 
we are interested in do not contain any phonon). The 
partial phonon matrix elements are calculated out, but 
we keep the general form for the magnetic states at this 
stage. In other words the phonons are integrated out 
and we end up with an effective matrix element acting 
between different magnetic states: 



(/|E.Pph|0) 

7 
6 



(/'|^7S..S,;+a + ^.(S, xS,+a) K7) 

ia 



E 



d,(D,.E) 



(8) 
(9) 



where Dg = '^j^qd^dsq=o is the amplitude of the instan- 
taneous electric dipole of the unit cell due to the phonon 
mode s with energy — i^g=o,s- The final magnetic 
state has an energy ujf. Qs = '^d add^'ds ^^'^ ampli- 
tude of the variation of the magnetic exchange energy due 
the atomic distortions of the phonon s (A^^ is the ampli- 
tude of the motion of the atom d, in the direction a due 
to the phonon s aX q — Q). Similarly, d" = X^/jd ^d'^-^ds is 
the amplitude of the instantaneous Dzyaloshinski-Moriya 
vector due to the phonon s. The resulting 7 and 8 depend 
on the bond considered. They would usually couple the 
nearest neighbors, but could be introduced for neighbors 
at larger distances if such super-exchange processes were 
likely to take place. They can be introduced on the basis 
of the symmetry which is usually reduced with respect 
to the crystal symmetry by the presence of the external 
electric field. Thus we have written an effective operator 
announced in eq. ^with A^-y = 
The selection rules are: 

• (i) D^.E 7^ 0: the virtual phonon s creates dis- 
tortions that carry an instantaneous electric dipole 
Ds. In other words, the phonon s must be infra-red 
active. 

• (ii) 

— (7s 7^ 0: The distortion of the unit cell due 
to the phonon s modulates the magnetic ex- 
change between the spins. The transition at 
^Stot = is allowed. 

— ds 7^ 0: It implies that the distortion of the 
unit cell due to the phonon s must break in- 
stantaneously the symmetry by inversion at 
the middle of the bond; so that to allow an in- 
stantaneous Dzyaloshinski-Moriya interaction 
which amplitude is given by ds. The transi- 
tions between states that differ by the spin, 
^Stot = 1, are allowed and have an intensity 

Suppose that there is only one phonon mode s which 
gives a major contribution to the sum. In addition, we 
know that this active phonon mode will appear in the 
infrared spectrum at the energy f2q=o,s, with an intensity 
given by Iph,s = (Ds-E)^. We can therefore rewrite the 
intensity of the A5tot = 1 line as: 



1 f^sds 



H 2 



2-/ 



iph,s 



(10) 



We denote by E the order of magnitude of the variation 
of the magnetic exchange energy due to the phonon and 
following MoriyajSiwe estimate dg ^ {^)E. That gives: 
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'ph,s 



(11) 
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This expression gives the intensity of such a process com- 
pared to the intensity of the optically active phonon. It 
is reduced by two factors: the spin-orbit coupling (in the 
cuprate materials, Ag/g can be 0.1) and the ratio of the 
energy modulation of the magnetic exchange due to the 
phonon by roughly the energy of the same phonon. The 
latter is difhcult to estimate : in CuGeOs, the first op- 
tical phonons have ^ lOmeV, and the modulation can 
be as large a.s E ^ ImeVm^ That gives /g ~ 10"**/^^. 

Another way to compare with is to consider that sin- 
glet excited states, as for example the S' = bound-state 
below the continuum in CuGeOa, appear in the optical 
spectrum due to the isotropic spin-phonon coupling (the 
7 term) . We denote their intensity by /|™9''=* . Then we 



have: 



jsmgiet ^ nieaus that if the sin- 



glet bound-state appears in the optical spectrum with an 
intensity J^^"9iet ^^le isotropic spin-phonon cou- 

pling, the triplet states should also appear with an in- 
tensity which is roughly 100 times smaller, if Moriya's 
estimate applies. 

Effect of a magnetic field. We consider a basic triplet 
excitation here. A magnetic field lifts the degeneracy of 
the triplet into three branches. When H || 5 (|| z), S*^ is a 
good quantum number and the transition should satisfy 
AS^ = 0. Therefore, only the mode = could be ob- 
served and its intensity does not depend on the strength 
of the field. By contrast, when the magnetic field is per- 
pendicular to 5, the wave-function is a superposition of 
wave-functions with different 5'^: 



an anisotropy in the g factor, which are both first order 
in the spin-orbit coupling, so that in the most favourable 
case (when no lattice selection rule forbids it), the ma- 
trix element is of order ^ Ag/g at best. In the second 
model, we consider an isotropic magnetic model, but we 
add the anisotropic spin-phonon coupling that we have 
considered above. The intensities of the transitions of 
the two models are given by: 



h-i = |(/|.gAi_Bh.Sf 



- [g^^Bhr ( ^ 



I, = |(/|D.E|0)P- ''^•'^^ ' 



So that the ratio is: 
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\giiBH 



9 J Wf-^l 



n,E 



(14) 



(16) 



where E — cH, c is the speed of light. D is given by 



D ^ eX where A 



is the amplitude of the motion 



of the ion and e is its charge. 
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(17) 
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0' _ Lf. ^' 



HI 



2\H^\ 
H 



|1,-1> 
11,1' 



(12) 
(13) 



where the vector notation stands for 13,3"). The transi- 
tion is allowed to the states ^E**' with quantum numbers 
S"*" = ±1 and the mode which energy does not depend on 
the field has no intensity. The magnetic field dependence 
is therefore very different from what is expected for mag- 
netic dipole transitions 1^ This is basically because the 
electric field conserves the 3^ quantum number. As we 
have just seen, however, in transverse magnetic field, 3^ 
is no longer conserved and the magnetic field-dependent 
branches may appear in the optical spectrum. They do 
indeed appear in CuGeOs 

We now compare the intensities of the magnetic dipole 
transitions with those of the electric dipole transitions 
that we have made explicit here. To make such a com- 
parison, we consider the following two models that give 
intensity to the optical transitions. First a purely mag- 
netic model and magnetic dipole transitions. In order 
to have an intensity, we need to add a static magnetic 
anisotropy, such as a Dzyaloshinski-Moriya interaction or 



with Mcu = 63g/mol (Mat ^ lO'^Skg), fl = lOmeV, we 
find A - O.lA. gfiB = 120/xeV/T. We take uj = 5meV for 
the energy of the magnetic mode and g = 2meV for the 
spin-phonon coupling. This estimation has to be taken 
with a pinch of salt because of the the crude order of mag- 
nitude given above, but it shows that there is no partic- 
ular reason to not consider the electric dipole transition 
due to dynamical Dzyaloshinski-Moriya interaction. 



III. APPLICATION TO CuGeOg AND SrCu2(B03)2. 

We compare the selection rules derived above with the 
experimental observation in CuGeOa. Experimentally, 
the absorption has been observed in the configuration 
E _L c but an extinction has been reported for E || Cfii 
even in the presence of a magnetic field. 

We have a natural interpretation of this fact: when 
E II c, the only contributions to S come from the vir- 
tual phonons s that have || c, or, in other words, 
the virtual phonons involved are those which are opti- 
cally active in this configuration. The vector S is given 



by E.rffAf^ 



ds^q—O 



where A^^^ q are the displacements 



of the atoms, the same as those that appear at higher 
energy in the real phonon state s. In the configuration 
E II c, the atoms in the phonon state s roughly move 
along the c-axis. In a crystal with many atoms per unit 
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FIG. 1: Three examples of distortions of the CU2O2 cluster 
with associated Dzyaloshinski-Moriya vectors. On the left, 
the motion along the y-axis creates a Dzyaloshinski-Moriya in- 
teraction whose vector is along z (the mirror plane CU2O2 and 
the one perpendicular which has the 0(2) atoms); on the mid- 
dle, the atoms move out of the plane and the Dzyaloshinski- 
Moriya vector is along y (mirror plane xz passing through 
the bond Cu-Cu) ; on the right, these distortions break the 
inversion center at the middle of the Cu — Cu bond. However, 
the two perpendicular mirror planes xy and xz imply that 
the Dzyaloshinski-Moriya interaction actually vanishes. For 
CuGeOs, the a;-axis is the c-axis and the xy plane is the plane 
of the Cu02 chains. 



cell, this is not exactly true and the displacements will 
acquire other components (a full study of the phonons 
that have been theoretically predicted in Ref. 30 does 
not change the picture). Then, according to the figure 
n the dynamical Dzyaloshinski-Moriya interaction is for- 
bidden = because of the mirror plane containing the 
atoms and the mirror plane perpendicular to the pre- 
vious plane and containing the Cu atoms. Therefore, 
the intensity vanishes in this special configuration. In 
other configurations, however, there is no such symme- 
try arguments leading to a cancellation of the dynamical 
Dzyaloshinski-Moriya interaction, and an intensity is ex- 
pected in agreement with the experiment performed in 
CuGeOa. 

We now consider the electric dipole transitions in 
SrCu2 (603)2 in greater detail. The obvious advantage 
of this compound is that, neglecting anisotropics, it is 
described by the Shastry-SutherlancSi Hamiltonian that 
possesses an exactly known ground state as a product of 
local singletsi^ Optical transitions have been observed 
between this ground state and each of the zero-field three- 
split triplet states^-'^® (see Fig. [JJ that have been de- 
scribed previouslyii^^ The probability of a transition be- 
tween the ground state ^'o and an excited state / is given 
by: 



(/l5^7S..S,-f<5,,.(S, xS,)|*o) 



(18) 
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We have restricted the operator He to the nearest neigh- 
bor spins (nn) in order to find the largest effect. The first 



FIG. 2: Excitation spectrum in SrCu2 (603 )2 f or two direc- 
tions of the external magnetic field (from Ref. l3a) . Definitions 
of the intensities of the optical transitions are also given. 



part of it does not change the total spin but may gener- 
ate transitions to the first excited states if the system has 
some anisotropy. We have considered previously the exis- 
tence of a Dzyaloshinski-Moriya interaction whose vector 
is perpendicular to the planei^ We have shown that such 
first-order anisotropy does not give intensity within the 
assumption of magnetic dipole transitions. Here we start 
by considering the electric dipole transitions generated 
by the first part of the operator H18() and in presence 
of the static Dzyaloshinski-Moriya interaction. Using a 
symmetry argument we show that this part actually van- 
ishes, is a conserved quantity so that we only need 
to consider the matrix elements with a 5^ = final 
state. The symmetry by the mirror plane perpendicu- 
lar to the (ab) plane and passing through a dimer is a 
symmetry of the crystal. In this symmetry, the ground 
state and the operator J2nn l^i-^j ^^''^ both even. How- 
ever the triplet state — adiabatically connected to 
the local triplet at J' = (the next nearest neighbor 
exchange) is odd. Then the matrix element vanishes. 
Additional spin anisotropy of the Dzyaloshinski-Moriya 
symmetry with extra in-plane components is present be- 
cause of the small buckling of the crystal structure at low 
temperaturesi^^ However, this (together with possible ex- 
change anisotropics) would, in any case, respect the same 
mirror-plane symmetry. So the first term is not expected 
to give intensity, because of this special symmetry. 

In SrCu2 (603)2, the transitions have been studied us- 
ing polarized electromagnetic waves and exhibit very pe- 
culiar polarisation properties: in the configuration E | 
{ab), at zero field, only the state at 24.2cm~^ (i.e. the 
=0 state [the middle state] ) appears in the spectrum, 
but an external in-plane magnetic field gives intensity in 
the two other modes (upper and lower modes) Simi- 
larly when the magnetic field lies in the (ab) plane, only 
the upper state at 25.4cm~^ (i.e. the = ±1) appears 
at zero magnetic field while an in-plane magnetic field 
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allows observation of the middle state, but not the lower 

oneM. 

We now show that these observations are compati- 
ble with the dynamical Dzyaloshinski-Moriya interaction 
which leads to the second part of the effective operator 
To explain these results we need to find the partic- 
ular pattern of dynamical Dzyaloshinski-Moriya vectors 
and then the Sij . That crucially depends on the direction 
of the electric field of the wave, according to eq. I^. In 
the following, we will determine the Sij but we restrict 
them to nearest neighbor interactions. 

Configuration E(t) \\ (ab). Let us consider first the 
case of a wave- vector of the electromagnetic wave paral- 
lel to the c-axis, then the electric field lies in the ab plane. 
According to the first selection rule (i), only the virtual 
phonons which carry an electric dipole \\ (ab) may 
contribute to the sum @. We basically assume that the 
main displacements of the atoms in such a virtual phonon 
mode are confined into the (ab) plane. We make the as- 
sumption that the main components of A^s are parallel 
to the electric field, so that we should be able to find 
the main components of the Dzyaloshinski-Moriya vec- 
tors dij,s (eq. IHJ). To estimate them (and then dij), we 
fix the atoms d at the distorted positions Xds and we then 
apply the Moriya's rules which give the constraints on the 
Dzyaloshinski-Moriya vectors. In this case, the plane re- 
mains instantaneously an approximate mirror plane for 
the crystal structure. Subsequently, the instantaneous d- 
vector between the spins, generated by the distortions, 
should be perpendicular to this plane (parallel to the c- 
axis). The effective operator is therefore written: 



E\\(ab) 



^<5,^(S,xS,r+^<5f(S,xS,)tl9) 



nn, A 



nn,B 



where z is here again the c-axis. We have introduced 
two different Sf'^ to take into account the existence of 
two dimers per unit-cell. Taking the same would not 
change the argument. In the following we take the no- 
tation Si = [{6ff + {Sff]/2. The operator (O does 
not break the symmetry by rotation around the c-axis. 
A transition to the = ±1 when the external magnetic 
field is parallel to the c-axis is still forbidden. Only the 
= triplet mode (at the middle of the others^) is 
allowed to appear in the spectrum (this is in agreement 
with the general symmetry argument given above since 
the electric field breaks the symmetry by mirror plane). 
This is in agreement with the experimental result at zero- 
fieldii^ We further predict that a magnetic field parallel 
to the c-axis does not change the picture and gives no in- 
tensity in the other branches. We can give an estimation 
of the intensity assuming an approximate wave-function 
for the excited state that we take from the strong dimer- 
ization limit. In this approximation, the excitation with 
S"^ = is a purely local triplet on the dimer A or B. This 
gives an intensity: 



\{^fll\HE\^o)\' 



mH\\) - 



(21) 



We now consider the effect of a transverse magnetic field 
(H _L c) on the intensities. A transverse magnetic field 
splits the modes into three branches (figure |21 left). To 
evaluate the intensity of each branch, we first calculate 
the excited states in the approximation used above, tak- 
ing into account the static Dzyaloshinski-Moriya interac- 
tion which is responsible for the zero-field splitting. Note 
that the other in-plane components do not play any role 
in the triplet spectrum at g = 0^ so that only the per- 
pendicular component appear in the following. 

The eigenvalues are in fact twice degenerate. The 
eigenvectors are denoted by ^'^^g ' and ^'^^g with en- 
ergies We then calculate the matrix elements as 
a function of the transverse magnetic field: 



^E\\{ab)y^^> 

We find: 



1 



2 1 + 
4 1 + 



(23) 
(24) 



where h = gfj,BH±/2D is the transverse magnetic field 
in the units of the static Dzyaloshinski-Moriya interac- 
tion. A transverse field transfers intensity into the lower 
and upper modes. The two curves given by l'^{H±) and 
I^{H±) + Ii^{H±) are shown in figure |31 together with 
the experimental results of Ref. [l^ We have used the 
non-renormalized value of Z? = 0.09meV extracted from 
the energy spectrumS (all the calculations we performed 
here are in the limit J' /J 0, so that we use the value 
of D we would have extracted from such a calculation 
and not the renormalized value). Note that if we take 
Ie{H±) and Ie{H±) for instance, they cross at a given 
field Hj_ = 2V2D/{giJ.B) ~ 2.1T, which is in good agree- 
ment with the crossing of the fitted intensities in the 
original experimental article {H± = 2.3T)i^. This is 
most probably coincidental since we are using the wave- 
functions that are not renormalized by the interaction 
J'. 

Configuration Eit) \\ c. We consider the case of an 
electric field perpendicular to the plane E(t) || c. Let us 
suppose that the atoms move out of plane. According to 
the figure^ the dynamical Dzyaloshinski-Moriya interac- 
tion would be in plane and perpendicular to the Cu — Cu 
bond. The dimers are, however, perpendicular to one 
another. Therefore the dynamical Dzyaloshinski-Moriya 
vectors of adjacent dimers should be perpendicular as 
well. The effective electric operator is: 



7 




10 12 



FIG. 3: Configuration E{t) \\ (ab). H || (ab). Intensity of 
the optical transitions from the ground state to the first split 
triplet state in SrCu2 (603)2 in the electric dipole approxi- 
mation {S II c). (+) + ( — ) is the sum of the intensity of the 
upper and lower mode and (0) is the middle mode. The the- 
oretical curves are given by eqs. 12311 and (1241 . There are 
no fitting parameters except an overall amplitude. The ex- 
perimental data are from Ref. US Note the intensity of the 
(-l-)-l-(-) mode is here twice as large as shown in the original 
experimental paperi— 



nn,A 



nn,B 



where 6 (respectively 5') is perpendicular to the Cu-Cu 
bond of the dimers A (resp. B), so parallel to y (resp. 
x). Note that we take the same \S\ and |^'|. Strictly 
speaking there is no reason why they should be the same 
but taking into account the special direction of the field 
we can reasonably assume that the motions of the atoms 
which belong to adjacent dimers are similar at least for 
the low-energy phonons. Let us apply this operator on 
the ground state which is approximately a product of 
singlet states on the dimers (we thus neglect the effect 
the static Dzyaloshinski-Moriya interactions have on the 
ground state which would give small corrections to the 
result). 



has been found in such a polarised configuration,-^ so 
that we conclude that D > 0. Only a detailed super- 
exchange calculation of D would be able to infer it. The 
matrix elements giving the intensities are given by: 



/+l+i(iJ||) ^ |(vI'+i^=+^|i/^,||.|vI'o)P - S^A (28) 
/-|-i(iJ||) ^ |(*-f ="'|i/ij||c|*o)P = (29) 

In zero external magnetic field, the two final states are 
degenerate so that the total intensity of the optical tran- 
sitions is the sum of the two, i.e. S'^/2. In a magnetic 
field parallel to the c-axis (z-axis), the upper mode splits 
into two branches with equal intensity (5^/4. 

Furthermore, we calculate the intensities as a function 
of a transverse magnetic field. The excited states ^'^^g''' 



and * 



(T,0)/ 
9=0 



are twice degenerate, so we calculate: 



E\\c 



iHi_) ^ \{^[t'o^\HE\^o)\' + |(vI/^S-^"|i/^,|vI/o)|^30) 



(T,0)/| 



We find the following expressions for the intensity of the 
upper (-I-), lower (-) and middle (0) states: 



8 [i + /i2] [±\/rT7?- 1] 

4 l + /l2 



2 (31) 
(32) 



where h — g^BH±/2D. The corresponding curves are 
given in the figure ^ Note that the crossing between 
and occurs at giiBH± = A^/2D, therefore at a 
field two times larger than in the configuration E || (ab). 
The agreement with the experiment is very good since 
such a balance of the intensities has been observedji^ 
The lower mode does not actually appears in the spec- 
trum experimentally and this is compatible with the low 
intensity we found. If we take the non-renormalized value 
oi D = 0.09meV, the crossing of the intensities occur at 
H± = 4.6T which is in good agreement with the exper- 
imental value (~ 6T), as well as the overall behavior of 
the curves. 



2V2 V 



A.S'=+1 
9=0 



9=0 

B,S^ = -1 
9=0 



-.S^ = -l 
9=0 



(26) 
(27) 



Note that ^'^Ig ^^"^ ^9^0 both eigenstates 

of the Hamiltonian restricted to triplet states with the 
same energy J -I- 2D. Depending on the sign of D, there- 
fore, only the upper mode or the lower mode should ap- 
pear in the spectrum. Experimentally, the upper mode 



IV. CONCLUSIONS 

In this paper, we have considered optical transitions 
with emission of one magnetic excitation, AStot — 1- We 
give a mechanism in terms of phonon-assisted transitions 
in which a virtual phonon is involved. The selection rules 
of such processes were made explicit: in brief we need a 
coupling to an infrared active phonon that breaks, at 
least instantaneously, the symmetry of inversion between 
magnetically coupled ions. The intensity of such a pro- 
cess has been estimated and we argue that it should be 
larger than a magnetic-dipole transition, at least in sys- 
tems in which spin-phonon couplings are appreciable. It 
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FIG. 4: Configuration E(t) || c. H || (afe). Intensity of the op- 
tical transitions from the ground state to the first split triplet 
state in SrCu2 (603)2 in the electric dipole approximation 
(<5 _L c). (+), (— ) and (0) are, respectively, the upper, lower 
and middle mode. The theoretical curves are given by eqs. 

and 13211 . There are no fitting parameters except an 
overall amplitude. The experimental data are from Ref. IT^ 

provides an alternative to purely electronic transitions 
that are not allowed when an inversion center is present. 

We note that we have considered uniquely the con- 
sequences of phonon assisted optical transitions in the 
context of single-phonon experiments, i.e. ESR and ab- 
sorption. The same mechanism can lead to processes in 
Raman scattering allowing single magnon creation, with 
similar selection rules concerning centers of inversion in 



the lattice. The effective operators will have similar sym- 
metry but are not identical, involving the polarisations 
of both incoming and outgoing photons. Experimentally 
there are extra contributions linear in both spin opera- 
tors and spin-orbit couplings that arc not present in the 
single photon case. While for the spectroscopy of sin- 
gle magnons in the materials studied, Raman scattering 
should be useful, single photon experiments may permit 
more direct comparison with microscopic estimates of in- 
tensities. 

In the final section we have studied the two specific case 
of CuGeOa and SrCu2 (603)2 for which polarised exper- 
iments are available. We have shown that predictions 
of the phonon-assisted theory agrees well both with ob- 
served extinctions and also, for the case of SrCu2 (603)2 
where detailed results are available, with the dependence 
of intensities as function of the external magnetic field. 
Further optical data should be analysed in terms of an 
effective operator of the Dzyaloshinski-Moriya symmetry 
for the matrix elements in the electric dipole approxima- 
tion. Potentially such optical experiments can provide 
a means of probing microscopically the spin-phonon cou- 
pling which may be relevant to other experiments, for ex- 
ample neutron inelastic scattering experiments at finite 
momentum transfer, and a way of studying four-spin cor- 
relation functions involving some sort of local chiralities. 
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